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Magnetization processes of spin-^ Heisenberg ladders are studied using strong-coupling expansions, 
numerical diagonalization of finite systems and a bosonization approach. We find that the mag- 
netization exhibits plateaux as a function of the applied field at certain rational fractions of the 
saturation value. Our main focus are ladders with 3 legs where plateaux with magnetization one 
third of the saturation value are shown to exist. 
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Recently there has been considerable interest in cou- 
pled Heisenberg antiferromagnetic (HAF) chains, so- 
called 'spin ladders', where one of the fascinating dis- 
coveries is that the appearance of a gap depends on the 
number of chains being even or odd (for recent reviews 
see e.g. Q). In this letter we study ladder systems at zero 
temperature in a strong uniform magnetic field. This is- 
sue has so far only been addressed for two coupled chains 
with a magnetization experiment on Cu2(C5Hi2N2)2Cl4 
H and theoretically using numerical diagonalization 
series expansions g and a bosonization approach 
These studies found a plateau at zero magnetization 
whose width is given by the spin gap in the otherwise 
smooth magnetization curve. In this letter we extend the 
theoretical approaches to three and more coupled chains 
using strong coupling expansions, numerical diagonaliza- 
tion and a bosonization approach. We find that in general 
the magnetization curves exhibit plateaux also at certain 
non-zero quantized values of the magnetization. 

To be precise, we concentrate on the zero-temperature 
behaviour of the following HAF spin ladder with N legs 
(kept fixed) and length L (taken to infinity): 

L N L 

i x—1 i—1 x—1 

-hJ2 S tx, (1) 

i,x 

where the Si ;X are spin-i operators and h is a dimension- 
less magnetic field. We assume periodic boundary condi- 
tions along the chains but investigate both open (OBC) 
and periodic boundary conditions (PBC) along the rungs. 
The magnetization (M) is given by the expectation value 
of the operator M = ^ Y,i, x S i,x- 

Our main result is that (M) as a function of h has 
plateaux at the quantized values 



^(i-(M))eZ. 



This condition also appears in the Lieb-Schultz-Mattis 
theorem Q and its generalizations [0J|]. There, it is re- 
lated to gapless non-magnetic excitations, but plateaux 
in magnetization curves appear if there is a gap to mag- 
netic excitations. 

However, there are other arguments which lead to (||) 
as the quantization condition for (M) at a plateau. A 
particularly simple one is given by the limit of strong 
coupling along the rungs J' 3> J (which has also proven 
useful in other respects §). At J = one has to deal 
only with Heisenberg chains of length N and the only 
possible values for the magnetization are precisely the 
solutions of (|): (M) € {-1, -1 + 2/N, . . . , 1 - 2/N, 1}. 

For N = 2 and J' 3> J this consideration predicts a 
plateau at (M) = 0. The boundary of this plateau is 
related to the spin gap simply by h c = A. The strong- 
coupling series for this gap reads fllOfl 
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where we have extended the result of Jl(| to fourth order 
(this further order is crucial to obtain a zero of the gap 
for J' > 0). 

For N — 3 and small magnetic fields one has a degener- 
acy which makes already first-order perturbation theory 
non-trivial. For OBC and \(M)\ < 1/3, the low-lying 
spectrum is then given by a spin-i chain in a magnetic 
field whose magnetization curve is well-studied (see pd|] 
and references therein). On the other hand, for PBC and 
\(M)\ < 1/3, the first-order low-energy effective Hamil- 
tonian for (|l|) at N = 3 turns out to be (see also Jl2| ) : 
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where the S x are su(2) operators acting in the spin-space 
and act on another two-dimensional space which 
comes from a degeneracy due to the permutational sym- 
metry of the chains. In particular, the usual spin-rotation 
symmetry is enlarged by an U(l) coming from the XY- 
type interaction in the space of the Fourier transforms 
along the rungs (first factor in (|J)). The Hamiltonian (ji|) 
encodes the magnetization curve in the strong-coupling 
limit of the three-leg ladder with PBC for \{M)\ < 1/3. 

At least for OBC, the situation is a little more 
favourable for the (M) = 1/3 plateau that one expects 
for N = 3. It turns out that one can use non-degenerate 
perturbation theory to compute the energy cost to flip 
a spin around this plateau and thus its lower and upper 
boundaries: 
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Finally, it is straightforward to exactly compute the 
upper critical field h uc for the transition to a fully mag- 
netized state for PBC, even L and arbitrary N; 



2 ( J + J') N even, 

(1 - cos (Tr^i)) J' + 2J N odd. 



(7) 



Actually, the result (Q) at N = 3, h uc = 3/2J' + 2J also 
applies to OBC. 

We now proceed with considerations that follow closely 
classical work on single Heisenberg spin chains |[l| and 
the 2D triangular HAF [Q. We have numerically cal- 
culated the lowest eigenvalues as a function of the mag- 
netization, wave vectors and coupling constants on finite 
systems with up to a total of 24 sites. 




FIG. 1. Magnetization curve for N = 3 at J'/ J = 3 with 
OBC. The thin full lines are for L = 8, the long dashed lines 
for L = 6 and the short dashed lines for L — 4. The thick 
full line indicates the expected form in the thermodynamic 
limit. The two diamonds denote the series (|^) and ([]) for the 
boundaries of the plateau. 



Fig. |] illustrates our main results with a magnetization 
curve for N = 3 at J' /J = 3 with OBC. The thin lines 
denote curves at finite system size and clearly exhibit a 
plateau with (M) = 1/3. The strong-coupling approxi- 
mations (||) and (g) for the boundaries of this plateau are 
in good agreement with the finite-size data. Finite-size 
effects are also small for the midpoints of the steps in 
the magnetization curves |llf| . From this one obtains an 
extrapolation to infinite L which is shown by the full line 
in Fig. |. 

A more compact representation is given by 'magnetic 
phase diagrams' which show the projection of the con- 
ventional magnetization curves as in Fig. [I] onto the axis 
of the magnetic field. We illustrate this in Fig. || with the 
case of a two-leg ladder which has already been studied 
in ||. On a finite lattice the possible values of (M) are 
quantized; the values of the magnetic field h/J where a 
transition from one such value occurs to the neighbour- 
ing one are shown by lines as a function of J'/ J. Regions 
without such lines denote plateaux in the magnetization 
curve while regions where they are close to each other 
correspond to smooth transitions in the thermodynamic 
limit. 
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FIG. 2. Magnetic phase diagram of the ladder with N = 2 
legs. The thin full lines are for L — 12, the long dashed lines 
for L — 10 and the short dashed lines for L = 8. The thick 
full line shows the gap (^). 

The magnetization curves of |^,|| with J' /J — 5 are 
located somewhat beyond the right border of Fig. 0. Our 
figure contains additional information about the depen- 
dence on J' I J. Note that the numerically obtained val- 
ues to full magnetization do indeed match with the ver- 
sion of (Q) for the ladder: h uc = 2J + J'. For the bound- 
ary of the plateau at (M) = one observes that finite-size 
effects are not substantial for J' /J > 3/2, and in this re- 
gion one observes also good agreement with the approx- 
imation (||) for the transition value of the magnetic field 
h c — A. At weak coupling, field theoretic arguments |l4| ] 
predict the gap to open linearly, A ~ J'. This is compat- 
ible with Fig. |^, though due to the large finite-size effects 
in this region much longer chains are needed to confirm 
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this linear behaviour conclusively ]l5[ |. 
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FIG. 3. Magnetic phase diagram of N = 3 chains with 
OBC. The thin full lines are for L — 8, the long dashed lines 
for L — 6 and the short dashed lines for L — 4. The two thick 
full lines denote (Q) and (^|), respectively. 

Figs. H and [| show the magnetic phase diagrams for 
N = 3 with OBC and PBC, respectively (Fig. is a 
section of Fig. || at J'/ J = 3). Both figures clearly ex- 
hibit a plateau with (M) = 1/3 at least in the region 
J' > 2J. The strong-coupling series (||) and (|J) for the 
boundaries of this plateau are also shown in Fig. ^ and 
in the aforementioned region J' > 2J, where finite-size 
effects are small, one observes good agreement between 
the expansions and the numerical results. 
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FIG. 4. Same as Fig. |3| but with PBC. 

Figs. H and [| differ at least in their details. For ex- 
ample, one observes that at the upper boundary of the 
(M) = 1/3 plateau in Fig. || it becomes favourable to flip 
two spins at once rather than one for J' > 3 J which may 
be interpreted as one signal of the frustration introduced 
into the system. 

At strong coupling the excitations above the 1/3- 
plateau in Fig. [| are described by ^ with all spins 



aligned along the field (S* = ^). This is an XY-chain 
and therefore massless, providing us with an example of 
a plateau in the magnetization curve with gapless non- 
magnetic excitations above it. 

The strong-coupling expansions as well as the numer- 
ical results obtained so far clearly show the existence of 
plateaux for sufficiently strong coupling. To see if the 
plateau persists in the weak-coupling region J'« J we 
use abelian bosonization (see e.g. 17|). The compu- 
tation to be presented below is similar to the ones per- 
formed recently in ||,[l8| , so we refer the interested reader 
to these two references for more details. 

The starting point for this weak-coupling expansion is 
the observation that at J' = one has N decoupled spin- 
2 HAFs in a magnetic field whose low-energy properties 
are described by a c = 1 Gaussian conformal field theory 
|pj| . This theory is characterized by the radius of com- 
pactification R((M),A) which depends on the magneti- 
zation as well as the XXZ-anisotropy A which we have 
here introduced for convenience. At zero magnetic field 
one has 0: i? 2 (0, A) = ^ (l - A cos -1 A). In general, 
this radius can be computed using the Bethe-Ansatz so- 
lution jl9| for the Heisenberg chain. For a qualitative 
understanding it may be helpful to notice that the ra- 
dius of compactification R decreases if either A becomes 
smaller or if the magnetization (M) increases. 

Now the low-energy properties of the Hamiltonian (|l|) 
at small coupling are described by the following Tomon- 
aga Hamiltonian with interaction terms: 

r r N 

i 

+ J2\ A 2 : cos (±k F x + V4tt(^ + : (8) 

i ^ 

+A 3 : cos (y47r(4>i - 0i+i)J : 



+A4 : cos 



with Hi = -^dxtfii, and A^ ~ J'/ J. It is convenient to 
rescale the fields by the radius of compactification R in 
(|), i.e. <j), -> y/A^Rcpi. 

For N = 3 chains with PBC we now change vari- 
ables from the fields 0i, (f>%, 03 to -7= (0i 

2), 



"2, <P3 W ^/f Vrl <" ' Y3)i 

^7= (<f>i - <h)y (y1 +02-203). Following Q one can 
show that the perturbation terms with coefficients A3 or 
A4 in (^) give rise to a mass for the latter two fields. 

Let us now consider the remaining field 0diag := 
■^YliLi^i- Radiative corrections to (||) generate the 
interaction term : cos(20di ag /-R) : 01 • For N — 3 and 
at (M) = this operator is irrelevant (in the region of 
A close to 1), which confirms that three chains weakly 
coupled in a periodic manner are massless. 
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Now we address the question of the appearance of 
plateaux for (M) ^ 0. Since this requires a gap for the 
(magnetic) excitations, we expect such a plateau to occur 
if the remaining field </>diag acquires a mass. In fact, the 
additional interaction term 

: cosOdiag/i?) : (9) 

survives in the continuum limit on a plateau pi"fl . The 
operator (|J) can appear only if (||) holds, as can easily 
been inferred from translational invariance of the origi- 
nal Hamiltonian ([!]) and the fact that k F = ^(1 — (M))ir 
and a one-site translation of the lattice Hamiltonian 
(|l|) translates into the internal symmetry transformation 

0diag — > <^diag + 2NRk F - 

For N = 3, (§) requires that (M) = 1/3. If one now 
estimates the radius of compactification R(^,l) follow- 
ing e.g. |Q, one finds that at J' = the operator (||) is 
slightly irrelevant for A = 1. The dimension of this oper- 
ator decreases with J' implying that the 1/3-plateau in 
Fig. [| extends down into the region of small J'. It should 
be noted that the appearance of the plateau for a given 
small J' crucially depends on the value of A, explaining 
why the numerical evidence in Fig. ^ is not conclusive 
in this region. For simplicity we have concentrated on 
PBC. The case of OBC is qualitatively similar but more 
subtle in the details and will be discussed in pi] ]. 

In this letter we have shown that spin ladders exhibit 
plateaux in their magnetization curves when subjected 
to strong magnetic fields. We have mainly concentrated 
on the plateau with (M) — 1/3 in three coupled chains, 
but also other rational values can be obtained by vary- 
ing the number of chains N . Quantum fluctuations {i.e. 
the choice of the spin S — ^) do play an important 
role - the plateaux disappear in general when one inserts 
classical spins into (Q). Nevertheless, by an alog y to the 
two-dimensional triangular antiferromagnet [P2fl it seems 
likely that such plateaux would appear in first-order spin- 
wave theory. Also with Ising spins the behaviour of (|lj) is 
different: At zero temperature the magnetization changes 
only discontinuously between plateaux values (i.e. no 
smooth transitions occur), and a gap opens both for odd 
and even numbers of chains. 

It would be highly interesting to check experimentally 
whether the plateaux can indeed be observed, in particu- 
lar since nowadays materials with a given number of legs 
can be engineered. The first non-trivial check would be 
to look for the (M) — 1/3 plateau in a three- leg lad- 
der. Here, the material Sr 2 Cu305 [^3| comes to mind, 
which is however not very well suited for these purposes 
since the necessary order of magnetic fields is not acces- 
sible today due to its large coupling constants. However, 
there are at least two-leg ladder materials such as the 
conventional (VO) 2 P 2 7 || or Cu 2 (C 5 Hi2N 2 )2Cl4 |] 
with much weaker coupling constants. A three-leg ana- 
logue of such materials could provide a testing ground 
for our predictions, in particular if such a material can 
be found with J' > 2 J where we would expect a clearly 



visible plateau in the magnetization curve at sufficiently 
low temperatures (cf. Fig. [l]). 
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